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Abstract. The high time resolution and large area of the Rossi X-ray Timing Explorer have been essential in the detection 
and characterization of high-frequency quasi-periodic variability in the flux from neutron stars in low-mass X-ray binaries. An 
unknown phenomenon prior to RXTE, kilohertz quasi-periodic oscillations (QPOs) have now been detected from more than 
twenty systems. Their high frequencies (up to 1330 Hz) imply that they are generated close to the neutron star, where general 
relativistic effects are expected to play an important role. I summarize current models for the kilohertz QPO phenomenon. In 
particular, I show that there is a significant domain of agreement among the models that can be used to constrain neutron star 
structure and look for signatures of highly curved spacetime in the properties of the QPOs. 



INTRODUCTION 

Neutron stars are important laboratories for physics at 
high densities. Unlike the matter in relativistic heavy-ion 
colliders, the matter in the cores of neutron stars has a 
thermal energy that is much less than its rest-mass en- 
ergy. Various researchers have speculated whether neu- 
tron star cores contain primarily nucleons, or whether 
degrees of freedom such as hyperons, quark matter, or 
strange matter are prevalent (see Lattimer & Prakash 
2001 for a recent review of high-density equations of 
state). In addition, the strongly curved spacetime around 
neutron stars implies that we could observe predicted ef- 
fects of strong gravity, such as frame-dragging or signa- 
tures of an innermost stable circular orbit. 

The fast timing phenomena observed from accret- 
ing neutron stars with the Rossi X-ray Timing Explorer 
(RXTE) offer outstanding opportunities for us to probe 
the regimes of high density and strong gravity. In par- 
ticular, the kilohertz quasi-periodic brightness oscilla- 
tions (kHz QPOs) observed from more than twenty neu- 
tron stars in low-mass X-ray binaries are promising, be- 
cause their high frequencies imply an origin near the star, 
where general relativity must play a role. 

Here we discuss briefly some of the proposals for 
the origin of these QPOs, and their implications for 
strong gravity and dense matter. Although currently no 
first-principles magnetohydrodynamic simulations pro- 
duce sharp QPOs (a situation expected to change in the 
next few years as computers become faster and more ef- 
fects can be included), we show that current observa- 
tional constraints on models are significant enough to 
allow fairly confident inferences. In § 2 we show that 
general relativity must inevitably influence the QPO phe- 
nomenon, independent of any detailed models. In § 3 we 



describe the constraints on models that follow from the 
observations, and the constraints on stellar mass and ra- 
dius that follow from the viable options. We conclude 
in § 4 by discussing what discoveries about dense matter 
and strong gravity can be expected with a future ~ 10 m 2 
timing instrument. 



THE EFFECTS OF GENERAL 
RELATIVITY 

Some of the most exciting potential implications of 
QPOs involve the effects of general relativity. These in- 
clude signatures of unstable circular orbits, and general 
relativistic frame-dragging. However, if general relativ- 
ity does not influence QPOs then the interpretations and 
implications are less clear. 

From the observational standpoint, a potential link to 
nonrelativistic systems was made by Mauche (2002), 
following the work of Psaltis, Belloni, & van der Klis 
(1999). These authors have shown that particular pairs of 
QPOs from sources (sometimes selected from more than 
two QPOs in a given system) follow a trend that, on a log- 
log plot, appears to link neutron star systems with black 
hole and white dwarf systems. Since white dwarfs are not 
significantly relativistic, Mauche (2002) concludes that 
the phenomenon as a whole cannot involve general rela- 
tivity, and hence favors a model such as the one proposed 
by Titarchuk and collaborators, in which the QPOs arise 
from classical disk oscillations (e.g., Titarchuk & Os- 
herovich 1999; Osherovich & Titarchuk 1999; Titarchuk 
& Osherovich 2000; Titarchuk 2002, 2003). Abramow- 
icz and Kluzniak have also proposed a mechanism of 
nonlinear disk resonances to explain black hole sys- 



TABLE 1. Upper peak frequencies and 
frequency ratios for black hole QPO pairs. 
GRS 1915+105 may also have a pair at «168 Hz 
and ks113 Hz (e.g., Remillard & McClintock 
2003). 



Source 


Vupper (Hz) 


Vupper/ Vl ower 


GRS 1915+105* 


67±5 


1.63±0.13 


XTE 11550-564+ 


272±20 


1.48±0.24 


GRO J1655-40" 


450±20 


1.5±0.13 



* Strohmayer 2001b 

t Miller et al. 2001; Remillard et al. 2002 
** Strohmayer 2001a; Remillard et al. 2002 



tems (Abramowicz & Kluzniak 2001; Abramowicz et al. 
2003a), which they suggest could extend to neutron stars 
as well (Abramowicz et al. 2003b). 

There are still many unknowns about the brightness 
oscillations from disks in binary systems, and it could be 
that there are some underlying mechanisms in common 
between black holes, neutron stars, and white dwarfs. 
However, whatever the detailed mechanism is, basic 
physical considerations require that general relativity 
will have an impact near neutron stars or black holes. 
For example, consider the highest frequency bright- 
ness oscillation ever observed, Vqpo = 1330 Hz, from 
4U 0614+091 (van Straaten et al. 2000). The frequency 
of the innermost stable circular orbit is < 1500 Hz for 
M>1.6M & (see below). This means that all frequen- 
cies (e.g., orbital, vertical epicyclic, and radial epicyclic) 
are altered by the curved spacetime. In turn, oscillation 
modes of the disk are altered as well. Even radiative 
transfer is affected subtlely, by light deflection effects. 
Thus, regardless of the underlying mechanism, general 
relativity will alter the basic picture. 

One can also look to observations to find that there 
are important differences between the types of sources. 
Abramowicz & Kluzniak (2001) made the prescient sug- 
gestion that for black hole QPOs there would be small 
integer ratios between frequencies, before the first detec- 
tions of such a ratio (see Table 1). There may also be 
sources with small integer frequency ratios different than 
3:2 (e.g., the possible 4:1 ratio found in 4U 1630^+7 by 
Klein-Wolt, Homan, & van der Klis 2003). This insight 
may well prove to be a key to understanding this phe- 
nomenon. However, neutron star systems do not have a 
similar preference for any particular frequency ratio. For 
example, the most recent data for Sco X-l (kindly pro- 
vided by Mariano Mendez) are plotted in Figure 1. As 
can be seen, the 3:2 ratio that is prominent in black hole 
sources is not evident here (the apparent clustering found 
by Abramowicz et al. 2003b used older, less precise data 
for Sco X-l). 
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FIGURE 1. Ratios of upper to lower QPO peak frequencies 
for the neutron star LMXB Sco X-l (data kindly provided by 
Mariano Mendez). The location of a 3:2 ratio is indicated by 
the dashed line. Unlike in the case of black hole QPO pairs, 
there is no preference in neutron star sources for a particular 
ratio or set of ratios. Similar trends are evident in other neutron 
star sources. 



In addition, there is strong evidence that the spin fre- 
quency plays a role in generating at least one of the 
two strong kHz QPOs observed in neutron stars. This 
evidence comes from a comparison of the frequency 
separation Av between kHz QPOs with the spin fre- 
quency v S pi n , in cases where both are known. In Ta- 
ble 2 we compare the range of observed frequency sep- 
arations with the spin frequency as inferred from the 
persistent oscillations in SAX J 1808-3658 (Wijnands & 
van der Klis 1998; Chakrabarty & Morgan 1998) and 
from the frequency of burst brightness oscillations in the 
other sources. The identification of the frequency of burst 
brightness oscillations with the spin frequency (or close 
to it) is well-established by observations of burst oscilla- 
tions in SAX J1808-3658 (Chakrabarty et al. 2003) and 
XTE J1814-338 (Strohmayer et al. 2003) that are ex- 
tremely close to the frequency of the oscillations during 
persistent emission. From Table 2 we see that all sources 
have Av close to v sp i n or v sp i n /2, and indeed 4U 1702- 
429, SAX J1808-3658, and KS 1731-260 have Av con- 
sistent with exactly v sp in or v sp i n /2, within the uncertain- 
ties. As indicated by Table 2, the frequency difference 
does vary in several sources, requiring modification of 
basic spin modulation ideas (see, e.g., Lamb & Miller 
2001), but the importance of the spin is clear. 

The involvement of the spin, e.g., through radiation ef- 
fects or stellar magnetic fields, indicates that there are at 
least some processes affecting neutron star QPOs that are 
different from those that generate black hole QPOs. No 



TABLE 2. Spin frequency and frequency 
separation for neutron star LMXBs. 4U 1702, 
SAXJ1808, andKS 1731 have single measure- 
ments of Av, with uncertainties indicated. 



Source 


Vspin (Hz) 


Av (Hz) 


4U 1916-053* 


270 


290-348 


4U 1702-429 t 


329 


333±5 


4U 1728-34** 


363 


342-363 


SAX J 1808-365 8* 


401 


195±6 


KS 1731-260 5 


524 


260±10 


4U 1636-536" 11 


581 


250-323 


4U 1608-5211 


620 


225-313 



* Boirin et al. 2000; Galloway et al. 2001 

I Markwardt, Strohmayer & Swank 1999 
** Strohmayer et al. 1996 

* Wijnands et al. 2003; Chakrabarty et al. 2003 

§ Smith, Morgan, & Bradt 1997; Wijnands & van 
der Klis 1997 

II Di Salvo, Mendez, & van der Klis 2003; Jonker, 
Mendez, & van der Klis 2002; Mendez, van der Klis, 
& van Paradijs 1998; Wijnands et al. 1997 

II Berger et al. 1996; Mendez et al. 1998; Mendez et 
al. 1999; Yu et al. 1997 



such relation is evident for white dwarfs (Mauche 2002), 
suggesting again that even if there is a master underlying 
mechanism at work, there are important differences be- 
tween the different classes of objects. General relativity 
does play a role for neutron stars and black holes, hence 
full understanding of the QPOs is promising for our un- 
derstanding of strong gravity effects. 

As an aside, the demonstration that Av « V sp i n /2 in 
SAX J1808-3658 is an indication of how the superb data 
available from RXTE, combined with sophisticated anal- 
ysis (e.g., Chakrabarty et al. 2003, Wijnands et al. 2003), 
are still facilitating qualitative leaps in our understand- 
ing. Prior to the SAX J 1808 analysis, I and many other 
researchers argued that Av ~ V sp ; n in all cases. The new 
results have forced modifications of the original models 
(for a recent proposal, see the sonic point and spin res- 
onance beat frequency model of Lamb & Miller 2003), 
proving again the importance of high-quality timing data. 



IMPLICATIONS OF QPOS 

As discussed in, e.g., Lamb & Miller (2003), observa- 
tions of kHz QPOs in neutron star LMXBs give a number 
of clues to their physical origin. 

• The spin frequency is involved in producing the 
observed frequency differences. The difference can 
be close to v sp ; n or V sp i n /2. So far, sources with 
Vspin > 400 Hz always have Av ~ V sp ; n /2, whereas 
sources with v sp i n < 400 Hz always have Av « V sp i n 
(Muno et al. 2001). 



• This appears to be a single sideband phenomenon. 
That is, if the spin frequency modulates some 
other frequency, only one additional strong QPO 
is produced (additional QPOs have been found in 
4U 1608-52, 4U 1728-34, and 4U 1636-536 by 
Jonker, Mendez, & van der Klis 2000, but these are 
much weaker than the primary peaks). This restricts 
models significantly; for example, amplitude mod- 
ulation of one frequency by another produces two 
sidebands of equal strength. 

• An excellent candidate for the other frequency is 
the orbital frequency or something close to it. The 
requirements are that the frequency be in the right 
range, while also being able to change frequency by 
hundreds of Hertz (for a review of the observational 
properties see, e.g., van der Klis 2000 or J. Swank, 
these proceedings). The orbital frequency has these 
properties. If the orbital frequency is involved, then 
because one expects accretion to align the stellar 
spin with the sense of the accretion disk over a 
time short compared to the accretion lifetime, the 
orbital and spin directions are the same and hence 
the orbital frequency is expected to be close to the 
upper peak frequency. 

Detailed models need to identify a mechanism that 
produces the QPOs, selects a particular orbital radius 
among many, and allows this radius to change signifi- 
cantly (as indicated by the changing QPO frequencies). 
The current leading candidates include some variant of 
a beat frequency model (e.g., Lamb & Miller 2003), or 
possibly a resonance with the spin, modulating other fre- 
quencies (D. Psaltis, presented at the "Neutron Stars on 
Fire" conference, Princeton, NJ, 11-13 May 2003). In the 
former case, the upper peak frequency V uppe r is identi- 
fied with a frequency close to an orbital frequency v or b 
at some special radius (e.g., the sonic radius; see Miller, 
Lamb, & Psaltis 1998), and in the latter case it is iden- 
tified with a vertical epicyclic frequency of a nearly cir- 
cular orbit Vverticai (e.g., Abramowicz et al. 2003a). For 
constraints on neutron star structure these amount to the 
same thing because V ve rticai ~ V or b outside a neutron star 
(for a discussion see Lamb & Miller 2003). There are de- 
tails of the observations (e.g., the conditions under which 
Av « V sp i n /2 instead of V sp i n ) that are not obvious from 
first principles (for some ideas see Lamb & Miller 2003), 
but the general constraints on models suffice to constrain 
masses and radii as long as V up per ~ V or b at some radius. 

Titarchuk and colleagues (e.g., Titarchuk 2003) have 
suggested instead that it is the lower peak frequency 
Viower tnat i s c l° se to V or b, with consequently differ- 
ent implications. In their model, the upper peak fre- 
quency Vupper is instead close to the hybrid frequency 
( v iower +4v n \g) 1 / 2 , where the magnetospheric frequency 
V m g ~ v S pin. This is an interesting suggestion, but the re- 



cent high-precision measurements of SAX J 1808 present 
puzzles for this model. Chakrabarty et al. (2003) show 
that the spin frequency is 401 Hz, rather than half this 
value. Wijnands et al. (2003) find a pair of QPOs, at 
V lower = 499 Hz and V up per = 694 Hz; the hybrid model 
would predict v U pp el - = (499 2 + 4401] 2 ) 1 / 2 = 945 Hz, 
in conflict with the observations. Similarly, data for 
KS 1731-260 present difficulties. Wijnands & van der 



Klis (1997) find Vi 



898 Hz and v. 



upper 



1159 Hz. 



The burst oscillation frequency is 524 Hz (Smith et al. 
1997). If v S pin = 524 Hz, the hybrid model predicts 
Vupper = 1380 Hz. If instead v sp i n = 262 Hz, the hybrid 
model predicts V U pp er = 1040 Hz. Both appear not in ac- 
cord with the data. For this reason, we will concentrate 



on models in which v. 



upper 



Vorb- 



In such models, measurement of V U pp er constrains the 
mass and radius of a neutron star. If is the radius of 
a circular orbit of frequency v or b> men clearly R < Sorb- 
in addition, the high quality factors of the QPOs require 
that they be produced outside the region of unstable cir- 
cular orbits predicted by general relativity. For a nonro- 
tating star, for which the exterior spacetime is described 
by the Schwarzschild geometry, the radius of the inner- 
most stable circular orbit (ISCO) is Risco — 6GM/c 2 . 
The constraint R 0Y \, > R places a mass-dependent limit 
on the radius; for example, for a nonrotating star R < 
(GM/4n 2 v 2 lb ) 1 / 3 (Miller et al. 1998). The additional 
constraint R^ > ^isco places an absolute upper limit 
on the mass and hence on the radius. When one consid- 
ers frame-dragging effects, the upper limits on the mass 
and radius are (Miller et al. 1998) 



M < 2.2 M Q (1000 Hz/ v or b ) (1 + 0.75j) 
R < 19.5 km(1000Hz/v orb )(l+ 0.2;) . 



(1) 



Here j = cJ/GM 2 is a dimensionless spin parameter, 
where J is the stellar angular momentum. If in a partic- 
ular case one believes that the observed frequency is in 
fact the orbital frequency at the ISCO, then the mass is 
equal to the upper limit given in equation (1). 

The highest frequency QPO so far detected with 
confidence has a frequency Vqpo = 1330 Hz (van 
Straaten et al. 2000), which would imply M<1.8M 
and R < 15 km for a system with spin parameter j = 0.1. 
These constraints essentially rule out the hardest equa- 
tions of state proposed (see Figure 2). The existence of 
the ISCO means that the frequencies cannot be arbitrar- 
ily high. If the radius at which the QPOs are generated 
gets close to the ISCO, a variety of signatures are possi- 
ble, including flattening in the observed relation between 
frequency and countrate, or sharp drops in the amplitude 
or coherence of the QPO (see Miller et al. 1998). Such 
signatures would confirm the presence of unstable or- 
bits, a key prediction of strong-gravity general relativity, 
and allow a direct mass measurement. It is possible that 
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FIGURE 2. Constraints from orbital frequencies. The 
1330 Hz curve is for the highest kilohertz quasi-periodic os- 
cillation frequency yet measured (for 4U 0614+091, by van 
Straaten et al. 2000). This curve is for a nonrotating star; the 
constraint wedge would be enlarged slightly for a rotating star 
(see Miller et al. 1998). The solid lines are mass-radius curves 
for different representative high-density equations of state. The 
mass-radius curves are all for equilibrium nonrotating stars; 
note that rotation only affects these curves to second order 
and higher. Curves Nl and N2 are for nucleonic equations of 
state; Nl is relatively soft (Friedman & Pandharipande 1981), 
whereas N2 includes significant three-body repulsion (Wiringa, 
Fiks, & Fabrocini 1988). PC has a sharp change to a Bose- 
Einstein condensate of pions in the core when the mass reaches 
~ I.8M0 (Pandharipande & Smith 1975). Equations of state 
Nl, N2, and PC are not modern (i.e., not fitted to the most 
current nuclear scattering data), but are included for easy com- 
parison to previous work on equation of state constraints. Curve 
S is for a strange star equation of state (Zdunik 2000). Curve Q 
is a quark matter equation of state with a Gaussian form factor 
and a diquark condensate (kindly provided by David Blaschke 
and Hovik Gregorian). 



the system 4U 1820-30 has already shown such a signal 
(Zhang et al. 1998), but there are complications in the 
spectral behavior that make this uncertain (Mendez et al. 
1999). 



PROSPECTS WITH A 10 M 2 
INSTRUMENT 

As discussed by Swank (these proceedings), a trend ev- 
ident from RXTE observations is that as the mass ac- 
cretion rate increases (inferred, e.g., from the S a index; 
Mendez et al. 1999), the frequencies of QPOs increase 
and the amplitudes of QPOs decrease. For many sources 
observed with RXTE, the inferred mass accretion rate 
continues to increase after the QPOs become unobserv- 



able. This suggests that a larger area timing instrument 
would be able to detect higher frequencies. For example, 
for 4U 0614+091, projection of the amplitude versus ac- 
cretion rate trends suggests that a 10 m 2 instrument could 
detect frequencies of ~ 1500 Hz (M. van der Klis, per- 
sonal communication). Other atoll sources could yield 
even higher frequencies. 

An important threshold is reached at ~ 1500 Hz, 
because this is at or above the orbital frequency at 
the ISCO for realistic masses. The orbital frequency 
at radius r in a Kerr spacetime of spin parameter j is 
D. = M'/VO 372 + ;'M 3 / 2 ) (e.g., Shapiro & Teukolsky 
1983, equation 12.7.19), and to R lsC o = 6M[l - 

(2/3) 3 / 2 ;'] (see, e.g., Miller & Lamb 1996), so the orbital 
frequency is 



Vorbjsco « 2199 Hz(M /M) 



11 



2 x 3/2 ■ 



(2) 



For j = 0.1, v orb ,isco < 1500 Hz for M > 1.58 M©. It 
has long been suspected that neutron stars in LMXBs 
are more massive than the canonical 1.4 M© because of 
mass transfer, and direct evidence of this has arrived 
recently. Nice, Splaver, & Stairs (2003) report that the 
22 ms pulsar in the 0.26 day binary J0751+1807 (with a 
low-mass white dwarf companion, likely a remnant after 
substantial mass transfer) has M > 1.6M e at better than 
95% confidence, and the mass could be well above this. 
Measurements of QPOs above 1500 Hz therefore have 
excellent prospects for stronger constraints on masses 
and radii, and even for detection of signatures of the 
ISCO. A QPO frequency as high as 1800 Hz would 
be large enough to argue against all standard nucleonic 
or hybrid quark matter equations of state, leaving only 
strange stars (see Figure 3). 

A qualitative advantage of a large-area timing instru- 
ment compared with RXTE is that there are a number 
of sources with strong enough QPOs that they could be 
detected in less than a coherence time. For example, in 
some states, the kHz QPOs in Sco X-l could be detected 
within «4 ms with a 10 m 2 instrument (M. van der Klis, 
personal communication). Because current detections are 
averaged over many coherence times, it is not possible to 
determine whether, e.g., QPOs are present at all times or 
whether they are superpositions of more coherent pulses. 
A larger area instrument would help resolve this, and if 
indeed there are underlying highly coherent pulses this 
could lead to substantial additional insights. For exam- 
ple, if the QPOs are caused by the orbits of radiating 
clumps, then observation within a coherence time would 
lead to detection of Doppler shifts, which when com- 
bined with the observed frequency would allow a unique 
solution of both the gravitational mass of the neutron star 
and the radius of the orbit. 
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FIGURE 3. Constraints on mass and radius for hypothetical 
detections of a 1500 Hz QPO and a 1800 Hz QPO, if they are 
identified with an orbital frequency. At 1500 Hz one expects 
signatures of the ISCO to be present; a detection of 1800 Hz 
would present strong difficulties for standard nucleonic equa- 
tions of state. The equation of state curves are as in Figure 2. 
Rotational effects are not included in this figure. 



In conclusion, RXTE observations have not only re- 
vealed a previously unsuspected phenomenon, but have 
constrained models of kHz QPOs significantly. With cur- 
rent data, we are just short of expected signatures of the 
innermost stable circular orbit, a crucial predicted char- 
acteristic of strong gravity. A larger-area instrument is 
likely to push us over this important threshold, and also 
to allow novel new methods of analysis that can detect 
qualitatively new phenomena such as periodic Doppler 
shifts of orbiting clumps. 
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